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Abstract 

We propose an efficient procedure to obtain Green’s functions by combining the shifted conjugate 
orthogonal conjugate gradient (shifted COCG) method with the nonequilibrium Green’s function 
(NEGF) method based on a real-space finite-difference (RSFD) approach. The bottleneck of the 
computation in the NFGF scheme is matrix inversion of the Hamiltonian including the self-energy 
terms of electrodes to obtain perturbed Green’s function in the transition region. This procedure 
first computes unperturbed Green’s functions and calculates perturbed Green’s functions from the 
unperturbed ones using a mathematically strict relation. Since the matrices to be inverted to 
obtain the unperturbed Green’s functions are sparse, complex-symmetric and shifted for a given 
set of sampling energy points, we can use the shifted GOGG method, in which once the Green’s 
function for a reference energy point has been calculated, the Green’s functions for the other energy 
points can be obtained with a moderate computational cost. We calculate the transport properties 
of a G6 o@( 10,10) carbon nanotube (GNT) peapod suspended by (10,10)CNTs as an example of a 
large-scale transport calculation. The proposed scheme opens the possibility of performing large- 
scale RSFD-NFGF transport calculations using massively parallel computers without the loss of 
accuracy originating from the incompleteness of the localized basis set. 


1 



I. INTRODUCTION 


Since the discovery of quantized conductance in atomic-size contacts, there has been con¬ 
siderable interest in studying the electronic transport in nanoscale systems. So far, a number 
of calculation methods to investigate the electron transport properties through a nanostruc¬ 
ture embedded by semi-inhnite electrodes have been proposed such as the nonequilibrium 
Green’s function (NEGF) method [1], recursion transfer matrix method [2], Lippmann- 
Schwinger method [3] and overbridging boundary-matching (OEM) method [HE]. Among 
them, the NEGF approach is most commonly used in connection with tight-binding mod¬ 
els and hrst-principles methods for solving transport problems. In the NEGF formalism, 
localized basis sets consisting of either atomic orbitals or Gaussians are extensively uti¬ 
lized because it is straightforward to partition a whole system into a transition region and 
electrodes. However, the incompleteness of the basis sets sometimes degrades the accuracy 
[6] and the diffuse functions in the basis sets prevent us from achieving high-performance 
computing on massively parallel computers. 

On the other hand, the calculation method using the real-space hnite-difference (RSFD) 
approach umn does not suffer from these problems. However, the linear-equation solver 
used to obtain the perturbed Green’s functions of a transition region becomes computa¬ 
tionally intensive as the numbers of grid points and sampling energy points increase. More 
efficient linear-equation solvers are required to perform large-scale hrst-principles transport 
calculations by the NEGF approach using the RSFD scheme. 

In this study, we propose an efficient scheme to compute perturbed Green’s functions 
that combines the procedure developed in the OEM method and the shifted conjugate or¬ 
thogonal conjugate gradient (GOGG) algorithm [91110]. In the NEGF method, the matrix 
to be inverted includes the self-energy terms of electrodes, which are energy-dependent and 
destroy the shift invariance of the Krylov subspace with respect to the energy of electrons. 
The total computational cost of obtaining perturbed Green’s functions is proportional to 
the number of energy points. On the other hand, unperturbed Green’s functions are hrst 
calculated by the OEM method and then the effect of the electrodes is introduced in the 
wavefunction-matching procedure. Since the matrix to be inverted does not include non¬ 
linear energy dependent terms, the unperturbed Green’s functions are obtained by solving 
a set of shifted linear equations with respect to the energy. The shifted linear equations 
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can be solved by the shifted COCG method, in which the calcnlation is actnally carried ont 
only for a reference energy point and the solntions for the other shifted energy points are 
obtained with a moderate compntational cost owing to the shift invariance of the Krylov 
snbspace. Employing the mathematically strict relationship between pertnrbed and unper- 
tnrbed Green’s fnnctions developed in the OEM method m, one can greatly decrease the 
compntational cost of the linear-eqnation solver. 

The rest of this paper is organized as follows. In Sec. the NEGF method nsing the 
RSFD approach is briefly introdnced. In Sec. in we state the problem of obtaining per¬ 
tnrbed Green’s fnnctions and introdnce the compntational scheme combining the procednre 
developed in the OEM method and the shifted GOGG algorithm together with a nnmerical 


example. In Sec. W, to demonstrate the efficiency and applicability of the newly devel¬ 
oped method, we present a transport calcnlation of a G6o@(10,10) carbon nanotnbe (GNT) 
peapod connected to GNT electrodes. Finally, we snmmarize onr work in Sec. |V| 


II. COMPUTATIONAL SCHEME 


A. Hamiltonian of whole system 


We consider ballistic transport throngh a nanostrnctnre sandwiched between two semi- 
inhnite crystalline electrodes as shown in Fig. The periodic bonndary condition is applied 
in the x and y directions, while the semi-inhnite bonndary condition is imposed in the z 
direction. The system nnder investigation can be divided into three regions: the left electrode 
(L), the transition region (T) and the right electrode (R). Ey introdncing localized basis sets 
or real-space grids, we can decompose a whole Hamiltonian H into the following matrix form: 



Hl 

Blt 

0 

H = 

-®lr 

Ht 

Btr 


0 

pt 

^TR 

Hr 


where Hi, Ht and Hu are the Hamiltonian matrices for the left electrode, transition region 
and right electrode, respectively. Bit {Btr) is the conpling term between the transition 
region and the left (right) electrode. Eecanse of electronic screening effects, the conpling 
between the left and right electrodes is zero when the transition region is snfficiently large. 
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We next derive the Kohn-Sham (KS) Hamiltonian [12] for the RSFD scheme. Since the 
direction of current flow is assumed to be along the ^ axis, the KS Hamiltonian is composed 
of the block-matrix elements with a dimension of N^y{= x Ny), where and Ny are the 
numbers of grid points in the x and y directions, respectively. Wavefunction values on the 
x-y plane at the point z = Zk are expressed by the column vector \I/(^fc) with N^y components 
given by {'ip{xi,yj, Zk) : i = 1,..., j = 1,..., Ny}. To simplify the explanation, we here 
restrict ourselves to the case of the central hnite-difference approximation, i.e., iV = 1 in 
Eq. (1) of Ref. (TJ and local pseudopotentials. The extensions to higher-order hnite-difference 
approximations and the inclusion of nonlocal parts of pseudopotentials are given in Ref. 0 
The KS equation for a whole system is expressed as a relation between three adjacent terms 
along the z direction 


- Bl^{zk-i) + [E - A{zk)\^{zk) - B^^{zk+i) = 0 (fc = -oo,..., -1, 0,1,..., oo), (2) 


where E is the energy of an electron, A{zk) is an Wiy-dimensional block tridiagonal matrix 
including the potentials on the x-y plane ai z = Zk and the coefficients of the hnite-difference 
approximation, and Bz is a constant matrix that is proportional to the iVa;y-dimensional unit 
matrix J. The discretized KS equation of the whole system is written as 


[E-H] 


^(^o) 

T(2;i) 




0 , 


( 3 ) 
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where H is the block tridiagonal matrix 



( 4 ) 


with Azf. = A{zk)- The partitioning in Eq. Q corresponds to those in Fig. [^and Eq. ([^. 
Since the KS Hamiltonian matrix for the RSFD scheme is very sparse and block tridiagonal, 
we can employ efficient iterative algorithms such as the steepest descent (SD) or conjugate 
gradient (CG) method to solve the eigenvalue and matrix inversion problems. 


B. NEGF method using RSFD approach 


This subsection gives a brief explanation of the NEGF method for later convenience. The 
NEGF method gives the electron transport properties without calculating the scattering 
wavefunction explicitly. The Green’s function of a whole system can be dehned as the 
resolvent of the Hamiltonian matrix of Eq. Q, 


G{Z) 




( 5 ) 


where Z{= E + irj) is a complex energy with r] being a positive inhnitesimal. The important 
quantity determining electron transport properties is the perturbed Green’s function of the 
transition region written as 


Gt{Z) 


Z-Ht- tdZ) - ±r{Z) , 


( 6 ) 


where 

= B[^gL{Z)BLT and Er(Z) = (7) 
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are the self-energy terms of the left and right electrodes, respectively. Here 


Gl^Z) = [Z- Hl]-^ and = [Z - Hr]-^ (8) 


are the surface Green’s functions of the semi-inhnite left and right electrodes, respectively. 
In the RSFD-NEGF scheme, since Brt {Brr) has only one nonzero iVa-y-dimensional block- 
matrix element B^ [see Eq. (|^], the self-energy terms are found to take the very simple 
forms of 


and 


where 


Sl(^) 


T.l{zo] Z) 0 ■■■ 0 

0 0 ■■■ 0 


0 0 ■■■ 0 


0 ■■■ 0 0 




0 ■■■ 0 0 

0 ■■■ 0 Eni^m+uZ) 


(9) 


( 10 ) 


Si(zo; Z) = BlgL{z-i,z_i; Z)B^ 

and 


( 11 ) 




( 12 ) 


with g{L,R}{zkT zi] Z) being the A^ 3 ,y-dimensional {k,l) block-matrix element of .ri( 2)- 
The perturbed Green’s function of the transition region Gt{Z) can be obtained by inverting 
the hnite-size matrix given as 


D{zo]Z)-'E l{zo]Z) -B^ 0 

-Bl D{zi-Z) -B, 

-Bl D{z^;Z) -B^ 

0 -Bl D{Zm+U Z) - 'LR{Zm+l, Z) 


( 13 ) 


where D{zi] Z) = Z — A{zi). 
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The expression for the conductance at a real energy E is given by the Fisher-Lee formula 


T{E) — Tr VE)G^{Zo, Zm+l] E)Tji{rn+l] E)G^{Zm+l, Zq] E) 
where G^{E) is the retarded Green’s function of the transition region defined as 


(14) 


G^j,{zk, zi; E) = \im GT{zk,zi;E + i7]), (15) 

G?-(£)(= |GJ(i!)l') is the advanced Green’s function, and ri(R) is the coupling matrix 
calculated from the imaginary part of the self-energy term, 

Tx{zk;E) = i[Ex{zk;E)-El^{zk;E)] {XeL,R). (16) 


III. EFFICIENT PROCEDURE TO OBTAIN PERTURBED GREEN’S FUNC¬ 
TION 


A. Difficulties in obtaining perturbed Green’s function by iterative method 


The simplest way to obtain Gt{Z) is the direct inversion of the matrix given by Eq. 


(13). However, direct matrix inversion becomes impractical when the matrix size increases. 


In addition, iterative solvers for linear equations, such as the SD and GG methods, are 


not efficient in this case because the matrix to be inverted given by Eq. (13) is neither 


Hermitian nor very sparse owing to the self-energy terms. In addition, the a computation of 
Gt{Z), which is carried out independently at each energy, is a computationally demanding 
task. To circumvent these difficulties, we propose a computational procedure employing 
the unperturbed Green’s function used in the OEM method in E] and the shifted GOGG 
algorithm PED] that uses the collinear residual theorems na in the following subsection. 


B. Unperturbed Green’s fnnction and shifted COCG method 

The unperturbed Green’s function of the transition is written as 

GriZ) = [Z - Ht]-\ (17) 
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The relationship between the unperturbed and perturbed Green’s functions is expressed as 
the following equation, which was derived in Ref. |TT1 


Gt{zo-, zi] Z) 


Z)Gt{zo, Zu Z) 

Gt{.Zi, zi; Z) 


0 



0 

Gt{zi, zi; Z) 

= Gt{Z) 

I 



0 

Grizni, zi; Z) 


0 

Grizm+i, Zi] Z) 


T^iZm+li Z^G’j'iyZYn+l-, Zi^ Z'j 


. the Ith 


(18) 


For I = 0 and m + 1, 


Grizo, zq] Z) — Qt{zo, zq] Z) I — Z)Qt{zo, zq; Z) 


-1 


GTi^ZYa+li Zq^ Z'^ I Z) 


-1 


X Grizm+i, Zo; Z) I-Y,^{zo] Z)Gt{zq,zo-, Z) 


^ -1 


Gt{Zq^ Zm+l] Z) — I — Qt{Zq^ Zq] Z)'^^{Zq] Z) Gt{Zo, Zm+l, Z) 


-1 


X 


I r( Zm+li Z^G'j'(yZin+li Zm+li Z^ 


-1 


(19) 


( 20 ) 


( 21 ) 


and 


1 -1 


GTi^ZYfi+li Z^ G'j'i^^ZYfi+li Z^^\^ Z^ I (Zfn+l 1 Z')Q'j’ (ZiYi^i , Z^n+x , Z') , (22) 

where Grizi, Zj; Z) is the modihed GT{.Zi^ Zj] Z) under the influence of the self-energy terms, 
ci-nd which is expressed as 

Gt{zo, zq] Z) = Gt{zq, zq; Z) 

~\~Gt{^Zq^ Zm+li Z) ^ ] p(Zfn-i-i^ Z^ 


X 


I GTi,Zm+l) Zm+1^ •^)5^ij('^»7i-l-lj Z') ^T('^m+1) Zqj Z'j 


-1 


(23) 


and 


^r(^m-i-i) ZjYi+ij Z'j Gt(^z^^i, Z'^ 

+GT{Zm+l, Zo', ■Z^)X)l(^o; Z) 


X 


I — Gt{Zo, Zq] -Z’)^^(^o; Z) Gt{Zo, Zjn+ 1 , Z). 


-1 


(24) 






















The advantage of computing the unperturbed Green’s function is that the matrix inverted 


in Eq. (17), Z — Hx, is very sparse and complex symmetric. Moreover, — Hx are 
shifted matrices for a given set of energy points Zj. {k = 1,2,...), which is in contrast to 
Zk — Hx — ’ZxiZk) — ^niZk) in Eq. (|^. Therefore, shifted methods such as the shifted 
COCG method [9] can be applied to obtain unperturbed Green’s functions. 

To solve the complex-symmetric linear equations {Zk — = Cj, where Cj denotes the 

jth unit vector, by the conventional GOGG method [12], we build up the Krylov subspaces 
independently at each complex energy Zk, which requires time-consuming matrix-vector 
operations for each energy. On the other hand, the shifted GOGG method makes use of the 
fact that Krylov subspaces are shift-invariant; to solve the linear equations for a given set of 
sampling energy points, the spanned Krylov subspace is common to all the energy points. 
Once the Krylov subspace has been built up for the reference energy point, we can obtain 
the Green’s functions for all the sampling energy points without matrix-vector operations 
during the iteration steps. 

It is noteworthy that we do not need to store the whole matrices of the unperturbed 


Green’s functions Gx{Z) in memory. To calculate the transport properties by Eq. (14), the 
block-matrix elements of the perturbed Green’s functions, Gx{zm+i, Zq] Z) and Gx{zo, Zm+i, Z), 
are required. In addition, the density of states is relevant to the diagonal elements of the 
perturbed Green’s functions, Gx{zi, zi; Z). The evaluation of Gx{zi, zi, Z) requires the ad¬ 
ditional computation of the off-diagonal block elements Gx{zo,zi,Z) and Gx{zm+i, zi, Z), 


which are analytically determined in an analogous manner to the derivation of Eqs. (19 )-(22) 


from Eq. (18). It is sufficient to store a couple of block-matrix elements of the unperturbed 
Green’s functions so as to perform matrix-matrix operations between the iVa-y-dimensional 
block-matrix elements of the Green’s functions and the self-energy terms. Thus, we can 
obtain the perturbed Green’s functions with a moderate computation time and memory 
load. 

There is no guarantee that the residual norms of all the sampling energy points satisfy the 
convergence criterion when that of the reference energy point is sufficiently small. Because 
the iterations are performed until all the residual norms satisfy the convergence criterion, the 
total number of iterations increases, resulting in the slight increase of GPU time in the shifted 
GOGG method with respect to the number of sampling energy points. In addition, when the 
norm of the residual vector of the seed system becomes small, the numerical precision of the 


9 







residual vectors for the sampling energy points degrades. The seed switching technique |T0] . 
which switches the reference energy point to the point where the residual norm is largest 
among the sampling energy points after the solution of the reference energy point converges, 
enables us to build up the residual vectors for the sampling energy points without loss of 
numerical precision and signihcant increase of the computational cost. The detail of the 
seed switching technique is introduced in Appendix. 


C. Numerical test 

To demonstrate the efficiency of the present scheme, we apply the shifted COCG method 
to the calculation of the unperturbed Green’s functions of the transition region for a Na 
atomic wire. The valence electron-ion interaction is described using norm-conserving pseu¬ 
dopotentials [in] generated by the scheme proposed by Troullier and Martins Exchange- 
correlation effects are treated within the local density approximation [TB] of density func¬ 
tional theory [iiiiin]. The central hnite-difference approximation is adopted for the second- 
order derivation arising from the kinetic energy operator in the KS equation. The dimensions 
of the supercell are 20.0 x 20.0 x 7.0 bohr and a single Na atom is contained in the supercell. 
The real-space grid spacing is chosen to be ~ 0.45 bohr, which corresponds to a total Hamil¬ 
tonian dimension of 12800. All calculations are carried out on a workstation with an Intel 
Xeon E5-2667v2 GPU. Figure]^ shows the GPU time versus the number of sampling energy 
points for the shifted GOGG and conventional GOGG methods. The initial reference energy 
point is chosen to be the Fermi energy. For the conventional GOGG method, the GPU time 
is proportional to the number of energy points because the Green’s functions are calculated 
independently at each energy point. On the other hand, in the shifted GOGG method, the 
matrix-vector operations are carried out only at the reference energy point and the cost 
of the scalar-vector operations to update the approximate solutions for all the shifted en¬ 
ergy points is negligibly small. The increase in the computational time is attributed to the 
time-consuming matrix-vector operations in the additional iteration steps to reach global 
convergence. 
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IV. APPLICATION 


We next compute the conductance spectrum of a C6o@(10,10)CNT peapod, in which a 
(10,10)CNT encapsulates a Cgo molecule, as an application of the proposed method to a 
large system. Figure shows the computational model. This system has been observed 
by transmission electron microscopy |2U] and its electronic structure has been intensively 
studied by theoretical calculations [2IH23]. However, the transport properties of peapods 
have mostly been investigated using the tight-binding approximation [23] because of the 
large system size required for calculations. 

The size of the transition region is taken to be 42.2 x 42.2 x 32.6 bohr and the transition 
region contains 340 carbon atoms. The initial carbon-carbon distance is set to 2.68 bohr 
[22] and structural relaxation is carried out. To determine the KS effective potential for 
the transport calculation, the conventional supercell indicated by the rectangle in Fig. 
is used, where the periodic boundary condition in imposed in all directions, and the grid 
spacing is ~ 0.45 bohr, which gives a matrix of 774144 x 774144 for the Hamiltonian of 
the transition region. The number of sampling energy points in the transport calculation is 
101. Calculations are implemented using 64 Intel Xeon E5-2697v2 CPUs. The total elapsed 
time required to calculate the unperturbed Green’s functions for 101 energy points is 15294 
second, in which the average number of iterations required to satisfy the convergence criterion 
is 6713, while the elapsed time and the average number of iterations are 6029 second and 
3122, respectively, when the unperturbed Green’s function is only calculated for the Fermi 
energy. 

The electronic band structures of the (10,10)GNT and G6o@(10,10)GNT peapod are 
shown in Figs. |^a) and l^b), respectively. In the band structure of the G6o@(10,10)GNT 
peapod, there are three flat bands above the Fermi level, which originate from the tiu orbitals 
of the isolated Ggo molecule. Owing to the hybridization between the orbitals of the Ggo 
molecule and tt orbitals of the (10,10)GNT, the vr bands of the (10,10)GNT are divided by 
the bands originating from the tiu orbitals. 

The conductance spectrum of the G6o@(10,10)GNT is plotted in Fig. Owing to the use 
of the shifted GOGG method, we can sample a large number of energy points so as to iden¬ 
tify the spiky dips in the spectrum. It is found that three dips appear near the energy range 
of the flat bands. Figure shows a real-space picture of the charge densities of scattering 
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wavefunctions, which is useful for understanding the spatial behaviors of the transport phe¬ 
nomenon. The charge density distribution with the energy at T{E) k, 2 spreads around the 
(10,10)CNT, while that with the energy at the dip of the conductance spectrum accumulates 
in the vicinity of the Cgo molecule. These results imply the incident electrons are scattered 
by the tiu orbitals of the encapsulated Ceo molecule. Since peapods have recently become 
experimentally accessible structures, this information should be helpful in understanding 
and predicting the resonant scattering in experimental transport measurements | 23 ]. 


V. CONCLUSION 

An efficient scheme to compute perturbed Green’s functions, which combines the proce¬ 
dure developed in the OBM method based on the RSFD approach and the shifted COCG al¬ 
gorithm, has been proposed. Since the number of grid points in the RSFD approach is larger 
than the number of bases in the tight-binding approach, the computational cost of obtaining 
the perturbed Green’s functions of the transition region has been the bottleneck. In this 
procedure, unperturbed Green’s functions are calculated and perturbed Green’s functions 
are obtained by employing the mathematically strict relationship between the perturbed and 
unperturbed Green’s functions. The main advantage of the present scheme for computing 
unperturbed Green’s functions is that the spanned Krylov subspace is common to all the 
sampling energy points. Therefore, once the Krylov subspace for the reference energy point 
has been built up, the Green’s functions for all the sampling energy points can be obtained 
without time-consuming matrix-vector operations. To demonstrate the potential power of 
the present method, the transport properties of a C6o@(10,10)CNT peapod suspended be¬ 
tween (10,10)CNTs are investigated. By sampling a large number of energy points owing to 
the use of the present method, the scattering attributed to the coupling between the states 
of the CNT and Geo can be observed in the calculated conductance spectrum. 

The present technique enables us to efficiently obtain perturbed Green’s functions in 
the framework of the RSFD approach. The RSFD scheme of first-principles calculations is 
free of problems concerning the completeness of the basis set that appear in the methods 
using localized basis sets of either atomic orbitals or Gaussians. Moreover, this scheme 
has the advantage of being scalable to massively parallel computers without compromising 
on precision. We believe that this development will be important for executing large-scale 
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transport calculations using massively parallel computers. 
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Appendix: Shifted COCG algorithm and seed switching technique 

The COCG method is a numerical method for complex-symmetric linear equations as 

Ax = b, (A.l) 

where A is not Hermitian but complex symmetric, A = ^ A^. Its algorithm gives the 

vectors x„, Pn and of the nth iteration as follows under the initial conditions, Xq = po = 0, 


: 1, and /3_i = 0 [TU]. 


Xfi Xji—i “f" iPn—1; 

(A.2) 

/-, . Pn—2^n—l A\ f^n—2^n—l 

= (1 + Q;„_iA)rn_i r„_ 2 , 

ttn-2 Oin-2 

(A.3) 

Pn Pn—lPn—1') 

(A.4) 

(r„i,r„i) 

(A.5) 

, _ {rn,rn) 

'n-l / N • 

(A.6) 


Now we solve m sets of the shifted linear equations (A -|- ail)x^’‘^ = b using the reference 
system (A -|- asl)x = b, where Uj {i = 1,2,... ,m) is a given set of the sampling energy 
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points and I is an identity matrix. It is proved that the residual vector for the sampling 
energy point and the residual vector for the reference system are collinear: 


= 


n 


Tn, 


where is a scalar variable [H]. Using the collinear relation of Eq. (A.7), the 
of the shifted COCG method is given as follows. 


(A.7) 

algorithm 


G.b 

n+l 

f^n—lCUn N (sA) /^n—lQ<n (s,i) 

cXyi—i (y.fi—1 

is,i) 

n 

<+i 

(A.8) 



(A.9) 

Pn-l 

(S,2) o 

- 

n 

(A.IO) 


(i) , (i) » 

= x\Bi + On-lPnll, 

(A.ll) 

Pn 

— ' n ^ Pn-lPn-li 

(A.12) 


where = 1 and the other initial conditions are same as the reference system. 

Using the information of the reference system, the shifted linear equations (A + = b 


are readily solved. Note that the computational costs for Eqs. (A.7)-(A.12) are negligible 
because they consist of scalar-scalar and/or scalar-vector products. In addition, and 
as well as the required elements of Xn'^ and pn^ for the sampling energy points are 
stored in memory during the iterations. The iterations are carried out until all norms of the 
residual vectors satisfy the convergence criterion. 

When the residual vector of the reference system is small, the numerical precision of 
the residual vectors of the sampling energy points degrades. The seed switching technique 
replaces the reference system by the energy point having the largest norm of the residual 
vector among the sampling energy points to circumvent this numerical difficulty. It is obvious 


from Eqs. (A.8)-(A.12) that and have to be updated after the reference system 


changes while an\ (3n \ Xn\ and pn'^ are independent from the choice of the reference system 
From the collinear relation of the residual vectors in Eq. (A.7), and 7i., 


are derived 


as 


= 


(s,i) 


{k = n-l and n) 

TTl ’ 


(A.13) 


14 














where s is the index of the new reference system. The extension to the switching to an 
arbitrary energy point cr(^ {cxi, (72,, <Jm}) is also available [TU] . 
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FIG. 1. Schematic of transport calculation model. The transition region is suspended between the 
left and right semi-infinite electrodes. 
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FIG. 2. CPU time required to obtain unperturbed Green’s functions of Na atomic wire. The red 
squares are the results obtained by the conventional COCG method and the black triangles are 
those obtained by the shifted COCG method. The energy points are shifted from the Fermi energy 
by 0.1 eV. 
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FIG. 3. Atomic structure of (10,10)CNT with single Cgo in transition region. The dashed lines are 
the boundaries of the supercell used to determine the atomic structure and KS effective potential. 
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FIG. 4. Electronic band structures of (a) (10,10)CNT and (b) C6o@(10,10)CNT. The Fermi level 
is marked by the dotted line. 
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FIG. 5. Conductance spectrum of C6o@(10,10)CNT. 
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FIG. 6. Charge density distributions of scattering wavefunctions of C6o@(10,10)CNT. (a) and (b) 
correspond to the energies indicated by the arrows in Fig. The spheres represent the positions 
of carbon atoms. Each contour represents twice or half the charge density of the adjacent contour 
lines. The lowest-density contour represents a density of 5.0 x 10“^ e/A^ 
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